Laser Cooling of Trapped Fermi Gases deeply below the Fermi Temperature 
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We study the collective Raman cooling of a polarized trapped Fermi gas in the Festina Lente regime, 
when the heating effects associated with photon reabsorptions are suppressed. We predict that by 
adjusting the spontaneous Raman emission rates and using appropriately designed anharmonic traps, 
temperatures of the order of 2.7% of the Fermi temperature can be achieved in 3D. 
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The realization of atomic Bose-Einstein condensation 
(BEC) [1] and degenerated Fermi gases [2-4], are among 
the most remarkable recent achievements of the Atomic 
and Condensed-Matter Physics. In the context of these 
experiments it has been recently proposed [5] that a well- 
known effect in superconductivity (see e.g. [6]), the BCS 
transition, could be observable in atomic Fermi gases. 
The experimental observation of such effect remains, 
however, as an open challenge. The BCS transition re- 
quires temperatures well below the Fermi temperature 
{Tp) and, unfortunately, present techniques, based on 
evaporative cooling in a Fermi gas with two constituents, 
have not yet been able to reach the required low tempera- 
tures. This is due to the Pauli blocking of the atom-atom 
collisions, which prolongs the cooling into time scales at 
which technical losses become dominant [7], Alterna- 
tive cooling mechanisms [8,9] combined with the exter- 
nal modification of the scattering length [10], have been 
proposed as possible mechanisms to achieve BCS. 

The aim of this Letter is to show that laser cooling 
offers a realistic and effective method to cool fermionic 
samples well below Tp. Evaporative cooling [2] involves 
the loss of a large fraction of the initial atoms, which in 
addition lowers Tp. On the contrary, laser cooling intro- 
duces in principle no atom losses, and therefore allows to 
cool larger fermionic samples, without significantly low- 
ering Tp during the process. In addition, laser cooling 
allows to cool polarized single-component Fermi gases, in 
which the atom-atom collisions are almost absent. Once 
such a polarized gas was sufficiently cooled, an interac- 
tion (for instance dipole-dipole one [11]) could be exter- 
nally induced, which can also lead to BCS transition [12]. 

Laser cooling of trapped bosons towards BEC have 
been predicted to work in the Festina Lente (FL) regime 
[13,14], in which the spontaneous emission 7 is smaller 
than the trap frequency u>. Experimental confirmation 
of this possibility has been recently reported in Ref. [15]. 
For fermions, apart from standard problems (reabsorp- 
tions) , laser cooling is obstacled by the inhibition of spon- 
taneous emission [16]. Such inhibition, combined with 
the already slow spontaneous emission rate in the FL 
regime, would lead to unacceptably long cooling times. 
This Letter shows that both, reabsorption problem and 



inhibition of the spontaneous emission, can be overcome 
in the Raman cooling of trapped Fermi gases. This is 
done first by adjusting the spontaneous Raman rate, and 
second by employing the anharmonicity of the trap. 

We consider N fermions with an accessible electronic 
three-level A scheme, with levels \g), e) and \r). We as- 
sume that the ground state \g) is coupled via a Raman 
transition to \e) (which is assumed metastable). Another 
laser couples |e) to the upper state r), which rapidly de- 
cays into \g). After the adiabatic elimination of |r), a 
two-level system is obtained, with an effective Rabi fre- 
quency f2, and an effective spontaneous emission rate 7. 
The latter can be controlled by modifying the coupling 
from |e) to \r). The atoms are confined in a non-isotropic 
dipole trap with frequencies u>!j. y z , ui* y z , different for 
the \g) and |e) states, and non-commensurable one with 
another. The latter assumption simplifies enormously 
the dynamics of the spontaneous emission processes in 
the FL limit. The cooling process consists of sequences 
of Raman pulses of appropriate frequencies, similar to 
those of Refs. [14]. We assume weak excitation, so that 
no significant population in |e) is present. This allows 
to adiabatically eliminate e), and consequently to con- 
sider only the density matrix p{t) describing all atoms in 
and being diagonal in the Fock representation corre- 
sponding to the bare trap levels. 

In principle, an inherent heating is introduced by 
the reabsorption of the spontaneously scattered pho- 
tons. Fortunately, it has been shown that such heating is 
largely reduced in the FL regime [13]. Due to the Pauli 
blocking of the s-wave scattering, the elastic collisions 
can be considered as absent, as well as the collisional two- 
body and three-body losses. Therefore, the main sources 
of losses are background collisions and photoassociation. 
For the considered laser intensities and atomic densities 
both loss mechanisms are negligible [17], in comparison 
to the losses introduced by the removal of excited-state 
atoms discussed below. 

In this Letter we use the notation of Refs. [14]. Us- 
ing the standard theory of quantum stochastic processes 
one can derive the quantum master equation which de- 
scribes the atom dynamics in the FL regime. One can 
then adiabatically eliminate the excited state, to obtain 
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the rate equations for the populations Nj in each level of 
the ground-state trap: 
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where the rates are of the form: 
•)2 r 27T 
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with 
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and A TO ; = J d(p J d6W(6, 4>)\r]i m (k L )\ 2 . In the above 
expressions 27 is the single-atom spontaneous emission 
rate, is the Rabi frequency associated with the atom 
transition and the laser field, rji m (ki,) = (e, l\e lkL ' r \g, m) 
are the Franck-Condon factors, W(9, (j>) is the fluores- 
cence dipole pattern, u-^ (wf) are the energies of the 
ground (excited) trap level to (I), and 5 is the laser detun- 
ing from the atomic transition [18]. We consider £1 < 7 
to ensure the validity of the adiabatic elimination of |e). 

Note that the rates (2) are nonlinear, due to their quan- 
tum statistical dependence on the number of atoms in 
each level. In particular, we clearly distinguish two differ- 
ent quantum-statistical contributions: (a) If the relevant 
levels are occupied, R m i vanishes, and the atoms remain 
for very long times in the excited state I, i.e. the spon- 
taneous emission is inhibited. This affects negatively the 
cooling process in two ways: first, the cooling times be- 
come very long, and second, the excited-ground collisions 
can occur and lead to heating and losses. In addition, the 
adiabatic elimination used in the derivation of Eq. (1) 
ceases to be valid, (b) The fermionic inhibition factor 
(1 — N n ) appears in the numerator of the rates, introduc- 
ing also a slowing-down of the cooling process. 

The above shortcomings can be overcome in the follow- 
ing way. First, note that 7 can be controlled at will in 
Raman cooling [13]. One can increase it gradually during 
the cooling in order to avoid the inhibition effects, but 
remaining in the FL regime. Still, even if one uses such 
approach some small fraction of the atoms will remain 
in the excited state after the cooling pulse, and has to 
be removed from the trap in order to avoid non-elastic 
collisions. The latter aim can be achieved by optically 
pumping the excited atoms to a third non-trapped level. 
This introduces a new loss mechanism which has to be 
taken into account in the simulations. 

In the following we shall consider a dipole trap charac- 
terized by a Lamb-Dicke parameter 77 = 2na/X — 2, with 
a = y/h/2muj being the size of the ground state of the 



trap, and A the laser wavelength. This is the case of cur- 
rent experiments in Mg [19], although our scheme could 
be also applied to more general traps, and other fermionic 
species. Let us first consider ID cooling, as a test of 
our ideas. Due to numerical demands in the calculation 
of the Franck-Condon factors for large quantum num- 
bers, we present the results for only N — 200 fermions, 
500 trap levels, and an initial thermal distribution with 
T = 0.65Tp, We stress, however, that the method could 
be employed for larger TV and an initial T > Tp. For Mg 
atoms and a laser wavelength A ~ 600nm, rj — 2 corre- 
sponds to a trap frequency w — 2tt x 5.7 kHz. A sequence 
of two pulses with respective dctunings S = —15, — 16a; 
was applied. When the system becomes degenerated, the 
atoms usually jump after a spontaneous emission from e) 
into levels of higher energy, since the lower energy levels 
are occupied. Therefore large negative detunings of the 
order of four recoil energies are required. Every cooling 
pulse was assumed to have fl = 0.25a; and a duration 
t p = lOOo; -1 . The spontaneous Raman rate was mod- 
ified during the cooling process in such a way that the 
averaged lifetime of the excited atoms was of the order 
of 0.05i p . This average was estimated by evaluating the 
rates (2). The atoms which did not return to \g) after 
the cooling pulse were removed from the trap. 

Fig. 1 (a) shows the time dependence of the tempera- 
ture, whereas in the inset the final averaged population 
distribution is depicted. Fig. 1 (b) shows the atomic 
losses due to removal of the excited atoms with inhib- 
ited spontaneous emission. As observed, T ~ O.lTp can 
be reached within 4 s. In order to achieve such low T 
with relatively small losses the value of 7 has been ad- 
justed from an initial value 0.5a; to 8a; in the final stages. 
Note that this does not imply a departure from the FL 
regime, since the inhibition effects drastically reduce the 
effective width of the excited levels. The temperature was 
calculated from the mean energy E with the expression 
E = 0.5NE F (l + (7r 2 /3) (T/7»2), valid for low T. 

To achieve even a lower T, the excitation of atoms 
deeply placed in the Fermi sea must be avoided. An 
experimentally feasible method to do it is to use an an- 
harmonic trap, in such a way that the Raman lasers are 
close to resonance just for the atoms within some de- 
sired band of energies; only those atoms are effectively 
excited. In particular, for strongly degenerated systems 
one should only excite a band of levels close to the Fermi 
surface. In the simulations we have used the dependence 
w^ 9 = wn(l — an) for the trap levels, where a is the 
anharmonicity parameter [20]. Sufficiently small anhar- 
monicity does not affect the Franck-Condon factors, but 
it does influence the resonance term in the denominator 
of the rates (2). The requirement that the pulse res- 
onant with the transition ui'f n+s — w^, affects only the 
atoms occupying the band of energies around to, put on 
a the following constrain: am max < l/4|s| with m max 
the maximum level considered in the simulation. This 
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inequality has to be fulfilled for each cooling pulse. The 
width S m of the affected band is related to 7. Compari- 
son of the two terms in the denominator of (2) leads to 
the following formula 7 ~ as8 m uj / (AR m ^ m+s ) [21]. 

We have considered a fcrmionic gas, already pre-cooled 
by the previous method, with an initial T = 0.087V and 
N = 181 (10% of the above considered initial number 200 
is lost during the first stage of cooling) . The trap param- 
eters were chosen as above, except for the addition of a 
small anharmonicity a = 1.04 x 10 . Such a, and a value 
7 = 0.02w were chosen to allow for transitions inside a 
band of about 20 levels. The cooling scheme consist of 
two pulses with detunings designed to be resonant with 
the transition from the Fermi level (181) and 47 and 48 
levels below respectively. Those large detunings are nec- 
essary to fill the holes deep inside the Fermi sea by the 
atoms from the surface. Both pulses have f2 = 0.016a; 
and a quite long duration t p = 5000a; -1 . The latter is 
necessary due to the slow Raman spontaneous process in 
a deeply degenerated system. The losses in this case are 
negligible (~ 1%) mainly due to the long pulses used. 

Fig. 2 (a) shows the time dependence of the tempera- 
ture, whereas its inset shows the final atom distribution. 
At the end of this cooling stage (after 5 s) the temper- 
ature decreases down to 0.017V. Fig. 2 (b) shows the 
atomic losses due to the removal of the excited atoms 
with inhibited spontaneous emission. 

Finally, we have analyzed the more interesting case 
of a 3D cooling. For simplicity we have considered an 
isotropic trap, in which the level energies have a de- 
pendence similar to the ID case [22]. We have taken 
N = 26235 atoms, 101 energy shells (176851 levels), 
and an initial thermal distribution with T = 1.147V. 
At T — the atoms fill 53 complete energy shells. In 
this case the method of increasing 7 during the du- 
ration of the pulses was not used, and thus a rather 
complex sequence of pulses was employed. The cool- 
ing was divided into 7 stages, each one of them con- 
sisting in cycles of two Raman pulses with detunings 
8/lu = {(-6,-7), (-14,-15), (-19,-20), (-19,-20), 
(-18,-19), (-20,-18), (-20,-18)}, and tuned to in- 
duce transitions from the bands centered at the shells 
{16, 40, 90, 70, 61, (53, 57), 56} with respective widths of 
{20, 28, 60, 40, 40, (40, 10), 20} levels. The spontaneous 
emission rates were chosen 7/0; = {(0.0075,0.0088), 
(0.025,0.026), (0.071,0.075), (0.048,0.05), (0.045,0.048), 
(0.48,0.028), (0.15,0.07)}. The detunings were chosen 
taking into account the 3D character of the trap states, 
since the excitation quantum number can only be de- 
creased in one direction in one cooling cycle. The first 
four pulses redistribute atoms occupying the levels below 
E F , making this part of the system more degenerate, next 
six pulses empty the trap levels well above E F and the 
last four affect the band around E F . We use f2 = O.87 
in each pulse, except in the 11-th where fi = O.I7, to 
ensure the validity of the adiabatic elimination. 



Fig. 3 shows: (a) the time dependence of T/T F during 
the cooling (solid line), and (b) the atom losses. The 
temperature was calculated with the help of the grand 
canonical expressions for the mean values of the particle 
number and energy. At the end of the process (f = 13s) 
the system reached T = 0.027TV. In the simulations we 
assumed that a fraction of the order of 1% of atoms is 
excited in every pulse, except sixth stage where there was 
about 7% excited atoms in every pulse. If such fraction 
is kept constant, similar cooling times are expected for 
larger N , the only limitation being the photoassociation 
losses. 

Our analysis can be extended to the case of a multi- 
component Fermi gas. In such a case, the different species 
can interact via s-wave scattering, although for a suffi- 
ciently degenerated gas such collisions should be largely 
suppressed, and in practice only possible close to the 
Fermi surface, inducing a perturbativc modification of 
the trap levels. The collisional dynamics can be split from 
the laser cooling dynamics [23] , and can be studied within 
the formalism of Quantum-Boltzmann-Master-Equation 
[23,24]. The collisions will introduce a thermalization ef- 
fect, which in fact should help in the cooling. The case 
of an unpolarized Fermi gas is particularly interesting 
when the interaction between different species is attrac- 
tive, since in such a case BCS pairing should be possible 
for temperatures below a critical value T c . For reason- 
able densities and scattering lengths, T c is well below Tp, 
and beyond the reach of the present evaporative cooling 
techniques. The methods discussed in our Letter offer an 
interesting perspective in this sense, since they allow to 
reach T <C T F . 

Polarized Fermi gases are by themselves also interest- 
ing, since once the gas had acquired a sufficiently low T, 
an interaction, such as for instance dipole-dipole one [11], 
could be externally induced. Such interaction, partially 
attractive, allows for the appearance of BCS transition 
[12]. BCS transition in a dipolar gas presents clear advan- 
tages compared to the one generated in a two-species gas, 
since for the latter the precise matching of both Fermi 
surfaces is essential, and experimentally very demanding. 

Finally, let us point that the full description of the 
laser-induced BCS transition within the Master Equa- 
tion formalism remains as a fascinating challenge. Such 
analysis will be the subject of future investigations. 
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FIG. 1. ID cooling by using the controlled Raman rate: 
(a) time dependence of the temperature and (b) losses. Inset: 
final averaged distribution. For details, see text. 
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FIG. 2. ID cooling by using the anharmonicity of the trap: 
(a) time dependence of the temperature and (b) losses. Inset: 
final averaged distribution. For details, see text. 
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FIG. 3. 3D cooling by using the anharmonicity of the 
trap: (a) time dependence of the temperature and (b) losses. 
Dashed lines separate different stages of the cooling process. 
Inset: final averaged distribution. For details, see text. 
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